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Introduction 

This workshop was sponsored by AIM and the NSF and it brought in participants from the US, 
Canada and the European Union to Palo Alto, CA to work to translate questions from quantum 
computing, complexity theory, statistical learning theory, signal processing, and data analysis to 
problems in geometry and representation theory. 

In all these areas, varieties in spaces of tensors that are invariant under linear changes of co- 
ordinates appear as central objects of study. Despite their different origins, there are striking 
similarities among the relevant varieties. We studied questions such as finding defining equations, 
hidden symmetries, and singularities. 

The AIM workshop was preceded by a graduate student workshop at MSRI on "Geometry and 
Representation Theory of Tensors for Computer Science, Statistics, and other areas," which was 
organized by Landsberg, Lim, and Morton. One goal of this workshop was to introduce the basic 
tools from geometry and representation theory which are used in studying problems that come up 
in these other areas. About half of the participants attended both workshops. 

At AIM, the mornings were devoted to introductory talks given by L. Valiant, P. Comon, G. 
Gour, J.-Y. Cai, L. De Lathauwer, and J. Morton. In the afternoons the work focused on the 
following 4 areas, however many participants worked with more than one group. 

(1) Multi-linear techniques in data analysis and signal processing 

(2) Geometric approaches to P?=NP 

(3) Matchgates and holographic algorithms 

(4) Entanglement in quantum information theory 

The participants split into these 4 groups, each of which contained a mix of practitioners and 
geometers. The groups then went to work on translating problems into a language which could 
be understood by both parties. After the translation work reached a stable point, the two goals 
were to determine what is already known in each area and to translate the open problems posed by 
the practitioners into well-posed problems in geometry. Finally, the participants made first steps 
toward answers - they determined simpler questions to start working on before the larger problems 
can be addressed; they made initial calculations to determine what should be true; and they made 
plans for collaboration focused on these problems. 
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Part 1. Signal Processing 

1. Report 

The decomposition of tensors into a sum of rank-one tensors has retained the attention of those 
in signal processing for quite a while. Collaboration with geometers continues to play an important 
role in determining what is already known and what tools can techniques may be immediately 
applicable to this and other problems in signal processing. 

An important aspect of the workshop was developing a dictionary to translate notions in signal 
processing to notions in geometry. This dictionary opened the lines of communication between the 
engineers and the geometers. We will recap some of the relevant definitions here. 

Let V\, . . . , V n be vector spaces (real or complex) and let V\®- ■ -®V n denote their tensor product. 
A tensor T is an element of V\ <8> • • • <8> V n . A tensor decomposition is the expression of a tensor 
as the linear combination of other tensors (presumably of lower rank) . The rank of a tensor is the 
minimum r such that T may be expressed as the sum of r rank-one tensors. 

The rank of a tensor is unchanged by scalar multiplication, so it is natural to work in projective 
space. The following are the geometric objects that we will need. 

The space of all rank-one tensors is the Segre variety, defined by the embedding, 

Seg : PVi x ■ ■ ■ x PV n -► P(Vi <g> ■ ■ ■ <g> V n ) 
([vi], . . . , [v n ]) i ^ [vi ® ■ ■ ■ ® v n ]. 

The space of all rank-one symmetric tensors is the Veronese variety, defined by the embedding, 

v d : PV -► P(S d V) 

H - lH d }- 

The variety of partially symmetric rank-one tensors is the Segre- Veronese variety defined by the 
embedding 

Seg du ..., dn :¥V 1 x---xVV n - P(S dl Fi ® • • • ® S d "V n ) 

(M,...,K]) ^ [( Wl ) di ®-..® («>„)**]. 

Suppose X C PV is a variety. The r th secant variety to X, denoted a r (X), is the Zariski closure 
of all embedded secant P r ^ 1 's to X, i.e., 



a r (X) = (J F(span{xi, . . . ,x r }) C PV, 

where the overline indicates Zariski closure. The subspace varieties, denoted Sub ri rn (Vi 
V n ) C P(Vi ® • • • ® V n ) are defined in [30] ch. 5. as 

'■<"» ® ■ ■ ■ » v„> - {pi * pcw ® ■ ■ ■ « v H ) i V ^f - C P( 1; " im( t ) r„ ) r " } 
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(Note: these are related to "rank varieties" in [58] ch. 7.) We also consider the symmetric subspace 
variety, 

Sub r (S n (V)) := {[T] G ¥(S n V) \ 3A C V, dim(A) = r, [T] C ¥(S n A)} C P(S n V). 

In this setting, we list various notions of rank. 

* The rank of the point [T] G P(Vi<g>- ■ -<8> V^) is the minimum r such that [T] can be expressed 
as the sum of r points on the Segre variety. 

* More generally, for an algebraic variety X C PV the X-rank of a point p G PV, denoted 
Rx{p) , is defined by Rx{p) '■= min{r | p G (xi, . . . ,x r ),x« G X}, where (.) denotes linear 
span. 

* The X -border rank is defined by R_x(p) '■= min{r | p G a r (X)}. In particular, a r (X) 
consists of all points with X-border rank no greater than r. 

* The symmetric rank of a tensor [T] G ¥(S n (V)) is the X-rank when X = v n (PV). 

* The partially symmetric rank of a tensor [T] € P(V / 1 ® <il <g> • • • ® V^® rfn ) is the X-rank when 
X = < Se< 7dl ,... A (Py 1 x ••• xPV n ). 

* The multilinear rank of a tensor [T] £ P(Vi <8> ■ ■ ■ <S> V n ) is the minimum (n, . . . , r n ) such 
that [T] G 5tt6 ri) ..., r „(Vi <g> • • • ® V n ). 

* The symmetric multilinear rank, (deonted X sym below) of a tensor [T] G P(5 ra y) is the 
minimum r such that [T] G Sub r (S n V). 

* When no X is specified, the convention is that the rank (border rank) of a tensor [T] G 
P(Vi (8) • • -®V n ) is the X-rank (respectively X-border rank) when X = Seg(¥V\ x • • • x PV^). 

* For X C PV, the generic X-rank is the X-rank of a generic tensor in PV, i.e. it is the 
smallest r such that a r (X) = PV. 

* A typical X-rank is an integer d such that the set {p G P(Vl (g) • • • <g) V n ) | Rx(p) = ^} has 
non-empty interior, i.e. it is an X-rank that occurs with nonzero probability. Note: over 
R there can be more than one typical X-rank. 

As far as applications are concerned, two important questions must be answered: (a) in what 
cases is a tensor decomposition unique, or at least when does there exist a finite number of such 
decompositions, and (b) how can one reasonably approximate a generic tensor by another of lower 
rank. A list of suggested problems from Comon (labeled P# below) were distributed before the 
workshop. During the workshop, the problems below were raised and discussed within the signal 
processing group. 

* The Alexander-Hirschowitz theorem and generic ranks: The Alexander-Hirschowitz 
theorem describes the dimension of ^(^(P^)), which is the expected one, with a list of 
exceptions. This theorem is important because, in particular, it determines the generic 
symmetric rank of a tensor (see [4-6]). G. Ottaviani gave a nice exposition of the proof of 
Alexander-Hirschowitz, (see also [21,46]). 

The extension of the Alexander-Hirschowitz theorem to not necessarily symmetric tensors 
should be the subject of future work (problem P5). It should be noted that progress in this 
direction has already been made, see [1]. Very recently, M. V. Catalisano, A.Geramita, and 
A.Gimigliano [12] have found that the secant varieties to Seg(F 1 x • • • x P 1 ) are all of the 
expected dimension with one exception. 

* Terracini's lemma: Terracini's lemma [52] is a tool which allows one to compute the 
generic rank numerically, and is especially useful when theoretical results do not allow one 
to obtain it. Terracini's lemma is also often used as a theoretical tool in the Alexander- 
Hirschowitz theorem, and the work of Chiantini-Ciliberto for example. 
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Terracini's lemma also holds true for tensors with Hermitian symmetries, which is useful 
in Signal Processing when the data are complex (e.g. in Telecommunications) and when 
the tensor is estimated from the data. An example was presented by P. Comon regarding 
cumulant tensors. In the complex case, and for a given order, there are several cumulant 
tensors that can be defined. An example of an order 4 cumulant tensor is 

Tijki = Cum{xi, Xj, x k , Xi }, 

where (*) indicates complex conjugation. This one has plain symmetries (ijkl — ► jilk) 
and Hermitian symmetries (ijkl — > klij or ijkl — > Ikji). Actually, these are the same 
symmetries as for the 4th order moment E{xi, Xj, x* k , x^}. See also [25]. 
Uniqueness of tensor decomposition: Uniqueness of a given tensor decomposition is 
related to papers of Mella [41], and Chiantini-Ciliberto, [16, 17]. In the cases where the 
Alexander-Hirschowitz theorem imply that there are at most a finite number of orbits, the 
work of Chiantini and Ciliberto gives a sufficient condition for uniqueness. Their condition 
is called weak defectiveness. This condition is straightforward, but to use it as a test is not 
a simple calculation. Mella gives an explicit list of some cases of uniqueness using different 
methods than Chiantini and Ciliberto. Certainly there is more work to do for the various 
types of tensor decompositions related to the other notions of rank mentioned above. 
Real tensors (P10): Tensor decomposition over R is more complicated. In particular, 
there can be more than one typical rank, but it is not known if there can be more than 
two typical ranks. This question has remained open after attempts by G. Ottaviani and P. 
Comon to find counterexamples. For the study of ranks in the real field, hyperdeterminants 
and other invariants should be useful. 

Non- negative tensors (Pll): Tensors with real positive coefficients are useful in spectral 
analysis and specifically to medical imaging. The approximation of a tensor by another of 
lower rank is well defined for tensors defined on the cone of non- negative reals R® n . A proof 
was sketched by L.-H. Lim. 

The approximation of a tensor by another of lower rank is also well-defined when the 
vectors in at least one mode are constrained to be orthogonal. A proof was sketched by L. 
De Lathauwer. 

A question that came up at the workshop is whether a specific metric is necessary. A 
metric is mandatory every time we have to deal with an approximation, which may happen 
for real positive, complex, or real tensors. However, it was claimed that its exact form is 
not necessarily essential. The question of decomposing a non-negative tensor into a sum 
of rank-one terms can be posed without a metric. Whether a metric would help to solve 
the problem is another question. It may be that a topology is sufficient, as for the complex 
case. 

Tensors with partial symmetries (P13): 

Consider the following types of tensors: 

(i) Tensors in the secant varieties to the Segre- Veronese, [T] € Cfc(<S'egi < i lj ... )< i Tl (Vi x • • • x 
V n )) C P(S dl Vi ® ■ • ■ ® S d "V n ). 

(ii) 3-tensors enjoying the circular symmetries, i.e after a choice of a basis {v l s } for 
each vector space V s we can write T = Tijkv\ (8) v 3 2 <8> S V\ ® V2 <8> V3 and the circular 
symmetries can be expressed as = Tj^i = T^ij. 

These tensors are encountered in particular when using the joint characteristic function, 
estimated from data measurements. These are also related to the BIOME and FOOBI 
algorithms. 
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(i) was discussed extensively during the workshop, whereas (ii) was only mentioned, but 
not focused on. It was pointed out that (i) can be useful for the decomposition of the 
characteristic function, see [23], as well as other types of partial symmetries. The example 
of circular symmetry quoted above was mentioned by L. DeLathauwer. See also [24] 
Kruskal's theorem for tensors: Given an explicit decomposition of a tensor T, Kruskal's 
theorem provides a test as to whether that decomposition is unique. A translation into al- 
gebraic geometry terminology given by J.M. Landsberg allows one to shorten the proof. 
Possible generalizations to symmetric and partially symmetric tensors as well as other ex- 
tensions were discussed. 

Comon's Conjecture (P15): Do rank and symmetric rank always coincide? In other 
words, do we have equality in the relation 

a k (» d (FV)) C a k (Seg(FV x • • • x PV))nP(S d V). 

This has been proved for at least two cases: (i) rank(T) < f-rankn h ,[k+i,...,d\(T) and 
k < d/2 (A proof of this fact, was derived and presented by D. Gross during the workshop. 
See Lemma 11.11 below for a definition of f-rank^ r .. i k],[k+l, ...,d](T) as well as the proof of 
Gross.), and (ii) under Kruskal's condition 2rank(T) < dn — d + 1, where n denotes the 
dimension. See also [21]. 

For other cases, the problem remains open, despite attempts by A. Bernardi and L. 
Oeding to find counterexamples. 

Rank versus border rank: We ask, for which varieties X is it true that the Rx(p) = 
R_x(p) for all p G P(V)? In general, rank < k is not a closed condition, so it is expected 
that rank and border rank rarely coincide, (cf. open problem 3). 

Note that for any variety X, the tangential variety t{X) is contained in the secant variety 
a(X). When X is the n-factor Segre variety, for instance, t(X) contains points of rank n, 
however all of these points have border rank 2. This example shows that there can be an 
arbitrary discrepancy between rank and border rank. 

Multilinear-ranks: L. DeLathauwer introduced a generalization of block decomposition 
of matrices to block decomposition of tensors. The group pointed out links between this 
and subspace varieties. The geometric objects that encapsulate this idea are the secant 
varieties to the subspace varieties. We ask, what is the generic X-rank when X is the 
subspace variety? This question comes down to computing the dimensions of the secant 
varieties a s (X). 

We also want to consider symmetric multilinear rank. We ask, what is the generic X- 
rank when X = Sub r {S k V)l Already when k = 3 these are important questions with great 
interest to engineers in signal processing. 

Strassen's (Generalized) Direct Sum Conjecture: Suppose X\ = Segi¥A\ x • • • x 
PA n ) and X 2 = Seg(PB 1 x • • • x FB n ), and let X = Seg(F(A 1 © B{) x • • • x F(A n © B n )). 
Suppose pi G A\ <S> ■ ■ ■ <8> A n , P2 G B\ <g) • • • <g) B n and let p = pi + pi- It is known that 
Rx{p) < Rxx{pi) + Rx 2 (P2), the conjecture is that equality holds. 

Note: The original conjecture was stated for bilinear maps, (i.e. the case n = 3), and 
Btirisser, et al. indicate that this conjecture is probably false (p360 [9]). In fact, Schonhage 
showed that the conjecture is false for border rank, so this seems to be a very subtle problem. 
P. Biirgisser gave a presentation of Schonhage's result during the workshop. 
Orbits (P6): For which vector spaces with group actions do there exist finitely many 
orbits? In the cases that there are finitely many orbits, what are those orbits explicitly? J. 
Weyman pointed out that this problem has been solved. A paper of V. Kac has addressed 
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the problem, [31], but there are mistakes, corrected in pg 523 [32], and in the multiplicity 
free case, [37]. 

• Explicit decompositions: Find explicit algorithms for the decomposition of tensors when 
the dimension is larger than two (c.f. (P7)). Several suboptimal algorithms exist, for 
example [3,24,25,35,36], however they are limited in that the symmetries are not fully 
exploited and the rank must be smaller than some upper bound, hence, these algorithms 
cannot decompose generic symmetric tensors. 

A suggestion was made by L. Gurvits to resort to doubly stochastic operators, but has 
not been pursued in depth. 

• Joint decompositions: Considering two tensors (of possibly different orders) simulta- 
neously may have two advantages: (i) There may be a possibility to restore uniqueness 
when the rank expected from the application exceeds the generic rank of the tensors (ii) 
There may be better robustness in the presence of noise. The problem (cf. (P18)) was 
discussed (J.M. Landsberg, P. Comon, J. Weyman, G. Ottaviani), but postponed to further 
collaborations. 

• Sylvester's theorem: This theorem was pointed out to be derived earlier by Gaspard de 
Prony (1755-1839), [26]. This is related to the signal processing problem of the estimation 
of frequencies of damped sinusoids. 

It should be mentioned that more generally, all the questions raised for Segre or Veronese varieties 
(genericity, uniqueness, etc.) can be posed for subspace varieties. 

For further reference for the ideas surrounding signal processing, tensor decompositions and 
questions of rank, we compiled the following list: [1,13-15,19,20,27,34,45,51]. 

1.1. From D. Gross. Let V = C n and Sym d (V) be the totally symmetric space in V® d . Let 

T G Sjm d (V). By a decomposition of T, we mean a set D = {a{\i, where af^ £ V for each i,k, 
and 

Oj = a\ 1] ® • • • ® af ] 

is a rank-one tensor that 

T = J>. 

i 

A decomposition D is symmetric, if all elements in D are symmetric. Let / C {1, . . . , d}. We write 

f-rankj jc (T) 

for the rank of T viewed as a degree-2 tensor obtained by grouping the respective indices in / and 
/ together (the 'f stands for "forgetful" or "flattening"). Also, if D is a decomposition, write 

^ (/) = {®aP' ) }* 

for the projection of the decomposition onto the factors in /. We get 

Lemma 1.1 (D. Gross). Let T e S d (V), d > 2. //, for any k, 

rankT < f-rank [1) ... )fc ] i[fc+1) ... )t q(T) 

then any minimal decomposition ofT is symmetric. 

What is more, let D = {a{\i be a decomposition of a symmetric tensor. Assume that for any 
I C {1, . . . ,d} where \I\ = d — 2 the set is linearly independent. Then D is symmetric. 
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Proof. We prove the second assertion first. Let / = {1, . . . , d — 2}. By assumption, is linearly 
independent, so there exists a dual basis {ctj}j fulfilling ai{aj) = 5ij for every a,j 6 -D (7) . Hence, 
for every i, 

a t (T) = a ! f- 1) ®af. 

Since ai{T) G Sym 2 (V), it follows that af ^ = af\ Now take I = 1, . . . , d — 3, d. Arguing as 

before, we arrive at af 2 ^ = af X \ By induction af^ = ap for all k,l, proving the claim. 

The first statement follows, because f-rank/^c (T) < rank(T) and equality trivially implies that 
is linearly independent for any minimal decomposition D. □ 

Note that the first assertion is formulated in terms of geometric properties of T, whereas the 
second statement just limits our ability to write down non-symmetric decompositions of symmetric 
tensors. 



2. Open Problems 

(1) Consider the case when X = Sub rj ... jr (V <g> ■ • • <8> V), and X sym = Sub r (S n (V)). Under what 
general conditions is it true that X-rank(T) = X sym -T&nk(T) for all tensors T in S n (V)? We 
can ask the same question for border rank. Notice that when r = 1, Subi i(V £3 ■ ■ -<S>V) = 
Seg(FV x • • • x FV) and Subi(S n (V)) = v n (V). Therefore, Comon's conjecture is a special 
case of this problem. 

It was remarked that for r = 1 some cases are known. For instance if dim{V) = 2 and 
arbitrary r, etc. see also, [22]. There are some very old results, for instance, Prony's method 
from 1795, which is equivalent to Sylvester's algorithm, [26]. 

(2) Consider X = Sub ri ^,^ rri (V\ ® • • • <g> V n ), and X sym = Sub r (S n (V)). What is the generic 
Jf-rank and X sym -iajok for X and X sym l Note that r±,...,r n , must have r; L < r^r^ for 

k distinct, otherwise it collapses. 
Suggestion of Weyman: One can relate this to a question of the existence of an open 
orbit. There is a close relation to the Castling Transform. Start with C ® C g <8> C r with 
p < q < r and r < pq. Then transform the question to C p <8) C q <8> C pq ~ r . Then you can 
restrict to an open set and get a correspondence on orbits. Therefore the property of the 
existence of an open orbit occurs simultaneously for (p, q, r) and (p, q,pq — r). 

(3) Let X = v k (FV) and consider Y r = {[T] € S k V \ X-rank(T) < r} with k > 2. Conjecture 
6.10 [21] claims that Y r ^ Y r for 1 < r < Rg, where Rg is the minimum r such that 
Y r = S k V. In particular, this says that the set of symmetric tensors of rank at most P is 
closed only for P = 1 and P = Rg. What can we say about this conjecture in this case and 
when X is a subspace variety or a Segre- Veronese variety? 

(4) What can we say about the uniqueness or finiteness of the corresponding decompositions 
t = t\ + £2 + • ■ ■ + td, where U G X, and X-rank(t) = d. We can ask this question when X 
is any of the varieties we have mentioned. Particularly interesting are the cases when X is 
a subspace variety, or a symmetric subspace variety. Again, many cases are already known, 
for instance see [51]. 

(5) Is the subspace variety version of Strassen's direct sum conjecture true? Let t\ G A\ <g) B\ <g> 
Ci, t 2 G A 2 ®B 2 <S)C2, so that h®t 2 G (A 1 © A 2 ) {B 1 B 2 ) ® (Ci © C 2 ) = i®B®C. Let 
Xj = Sub a r tb i tC r(Ai®Bi®Ci) C P(A<g)Bj(8)Ci) and J = Sub a ^ c {A®B® C) C P(A(8)5(8)C) 
be the respective subspace subvarieties. We ask if the following equality holds: 



X-rank(ti © t 2 ) = AVrank(ti) + A" 2 -rank(t 2 ). 



8 



LUKE OEDING 



What, if any, relationship must there be between X, X\ and X<{! What is the status of 
the conjecture for other varieties X,X\,XiP. This question can be posed for an arbitrary 
number of factors as well. 

(6) The previous 4 questions are posed over the base field C. Over R, what are the typical 
V-rank's and X SJ/m -rank's in the semi-algebraic sense? 

(7) Tensor product of cones over orthants: Let A + = R+, B + = R+,C+ = R+. Let X + C 
A+ ®M + B + ®k + C + =: V + be a semi-algebraic variety. Study the V + -rank. A tensor 
with all non-negative coefficients can be written as a non-negative linear combination of 
tensors of lower rank. In this case we must have X + = Seg+ = (<S'n6i ) i i i)+, and X + = 
Seg(¥A x Pi? x PC) n V+. It was remarked that (of course) this rank can be larger than 
the X-rank. 

Part 2. Mulmuley-Sohoni Approach to VP vs. VP vs. VNP 

3. Report 

The Mulmuley and Sohoni approach to P vs. NP has as its primary goal to "reduce hard nonexis- 
tence problems in complexity theory to tractable existence problems in geometry and representation 
theory" , [44] . Before we get to this geometry and representation theory, we will briefly review some 
of the relevant notions. 

During the workshop, L. Valiant described algebraic versions of the complexity classes P and NP 
(these are commonly referred to as VP and VNP), [56]. Computing the permanent is known to be 
an VNP-hard problem, however, by using techniques from linear algebra, the determinant can be 
computed in polynomial time. Valiant's approach to P vs. NP suggests attempting to show that 
if one computes the permanent of an m x m as the determinant of an n x n matrix, that n must 
grow faster than a polynomial in m. (For reference, see [42,43].) 

P. Biirgisser, J.-Y. Cai pointed out that Toda defined a skew arithmetic circuit (see [53-55]). 
An arithmetic circuit can be described by a class of finite directed acyclic graphs. All of the edges 
(usually called wires) are directed towards a single output node. Each node has either zero or two 
incoming wires. Nodes with no incoming wires are called inputs and the input to each wire is either 
a constant or a variable. Nodes that have two incoming wires are called gates, and each gate does 
a single elementary computations on those inputs. 

In an arithmetic circuit, there are only two types of gates allowed - multiplication and addition. 
So each gate takes precisely 2 polynomials as inputs and outputs either the sum or the product of 
the input polynomials. A gate may have many identical outputs, but there is only one gate that 
has no outgoing wire, and this contains the final output of the circuit. A skew arithmetic circuit 
carries the additional requirement that every multiplication gate have at least one of its inputs 
either a variable or a constant. Locally, a skew arithmetic circuit for multiplication looks like the 
following. 

Xi g 
\ / 
(*) 
I 

The inputs of this circuit are the variable Xi and the polynomial g and the output is their product, 
Xig. The final output of an arithmetic circuit is a polynomial. Toda showed that the determinant 
of an n x n matrix can be computed by a skew arithmetic circuit whose size is bounded by a 
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polynomial in n. A sequence of polynomials is in VP if it can be computed by a sequence of skew 
arithmetic circuits whose sizes and degrees are bounded by a polynomial in the number of inputs. 

Toda also defined a notion of weakly skew. Malod and Portier, [38] recall that a weakly skew 
arithmetic circuit is an arithmetic circuit with the following additional requirement: of the two 
sub-circuits above each multiplication gate, at least one must become disjoint from the rest of the 
circuit upon deletion of its wire to the multiplication gate. This may be depicted roughly as follows: 

□ — disjoint — □ 

\ ' / 
(*) 
I 

The boxes represent sub-circuits of the original circuit and (*) is a multiplication gate. A sequence 
of polynomials is in VPws if it can be computed by a sequence of weakly skew arithmetic circuits 
whose sizes and degrees are bounded by a polynomial in the number of inputs. These are examples 
of Valiant's complexity classes. It is not known whether VP = VPws i n general, however, Valiant 
comments that if you allow nlog(n) - i.e. quasi-polynomial, then all of these classes collapse. 

The sequence (det n ) is complete for the class VPws [38]. Since the sequence (perm m ) is VNP 
complete, Mulmuley and Sohoni are trying to prove that VPws 7^ VNP, where the class VPws 
is defined analogously to VP as in the paper of Biirgisser, [8]. 

Since we don't want to be deceived by bad choices of coordinates, we work with the appropriate 
orbits of perm and det. The Mulmuley-Sohoni approach to P vs NP is to determine if there exists 
a polynomial n(m) such that 

[l n ~ m perm m ] G [GL{W).det n ], 

where W = Mat nxn C and I is a linear form. We also want to know what we can say about the 
inclusion 

[End{W).det n ] ^ [GL(W).det n ]. 
It should be noted that the description of the orbit closure of a given polynomial is in general very 
complicated. 

The strategy suggested is to understand the homogeneous coordinate rings. We should try to un- 
derstand the modules of functions S C C[GL(W) .l n ~ m perm m ] which do not occur in C[GL(W).det n ]. 
Mulmuley-Sohoni call such an 5 a representation theoretic obstruction for l n ~ m perm m to lie in the 
variety GL(W).det n , [44]. 

Remark 3.1. Let V be a G- module for a reductive group G, and let v € PV. The we have 
the following inclusions G.v C G.v C PV. The big problem is to understand the restriction 
C[G.t>] <^-> C[(jt.i>]. The only general tool is a theorem of Kostant. It says that if the closure G.v is 
normal and G.v \ G.v has codimension 2, then the restriction map is an isomorphism. However in 
our examples these conditions are almost certainly not satisfied. Representation theory gives the 
following decomposition: Let G v denote the stabilizer in G of v, then 

C[G.v] = C[G/G V ] = V A * <S> V A G " . 

XeIrrep(G) 

This gives information about the coordinate ring of the orbit closure. 

Step 1: The first example is if v = [det n ], and W = M n C = A <g) B, then 

C[GL(W).det n ] = K x>S n, S nS x (A ® B), 
where K\^^ n are Kronecker coefficients. 
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The alternate description is that 

C[GL{W).det n ] = Sym{C p ®C n ® c n^SL n xSL n 

If A has 2 parts (i.e. p = 2), then the coefficients K\^ n ,S n are known. This is because the ring of 
invariants in this case is well known to be isomorphic to the ring of semi-invariants of a Kronecker 
quiver, see [50]. 

The connection is that C p ®C q <g>C r can be viewed as any of the following: (i) p copies of C q ®C r . 
(ii) p maps from C q to C r . (iii) the representation space of a quiver with 2 vertices 1,2 and p arrows 
from vertex 1 to vertex 2, of dimension vector (q,r). This is a generalized Kronecker quiver, for 2 
arrows we get just the Kronecker quiver. If we look at occurrences of S\C P <g> S^C 9 <g> S u C r in there, 
with A having 2 parts, the same representation occurs with the same multiplicity in the coordinate 
ring, as for p = 2. 

The Kronecker coefficients can be calculated using the Sylvester formula for the covariants of 
binary forms. Computing the Kronecker coefficients in general is a very difficult problem. For 
the computational complexity of these computations, see [10]. Immediately after the workshop, 
Manivel was able to make progress on this problem, [40], (cf. problem (4) below). 

There are special cases when, K^ ^, are known to be computable in polynomial time. They 
are as follows. (1) one of partitions has one row. (2) two of the partitions have two rows. (3) two 
partitions are hooks. (4) one partition is a hook and another has two rows. (For these results and 
more, see [47-49]). 

We want to find conditions which will determine when K Xj s n ,S n is zero or n °t- A central question 
is, "Is this decidability problem in P?" Mumuley conjectures that this decidability problem is in 
P. However, this is an open question, [42]. Another possibly interesting article on this topic is [18]. 

Computational efforts have determined that K\ t s n ,S n is almost never zero. In fact Derksen and 
Landsberg compared what occurred in the symmetric algebra in each degree. They did this in 
order to have a first approximation as to which modules could even be modules of polynomials. 
Unfortunately, not only did each module occur, but with high multiplicity. 

Ultimately, the success of this line of reasoning will depend on asymptotic questions. It becomes 
necessary to describe a certain polytope which is asymptotically nonzero. There is a hope that this 
method will be reasonable because the K\^n ^ n are not just generic Kronecker coefficients, they 
have a special format. 

The idea was advanced to use the moment map to establish the vertices and/or faces of the cone 

C{p,q,r) = {(\,n,u) | K Xtli , v ^0} 
This is an interesting problem, and L. Manivel described what is known during the workshop. 
Step 2: The second example is the permanent v = perm m . We need to consider 

C[GL(W).perm m ] = M X S X (A ® B), 

where M\ = Kx ifl)U , and we require that \i (resp. v) is such that [S^A]®" / ( resp. 

[S u B]f n ^ 0). This is seen to be a more tractable question, and there are a few papers on this. 
For n = 3, Weyman could show the formula for [S^AJ as an ©3-module. 

Step 3: The final question is what happens when we expand to more variables, i.e. when making 
the transition from perm m to l n ~ m perm m . A polynomial / is called stable under the action of 
some group G whenever G.f is a closed set in the Zariski topology. Unfortunately, while perm m 
is stable under the action of SL m 2, the polynomial l n ~ m perm m is not stable under the action of 
SL n 2 (for a short exposition of this concept, see [2]). 
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Mulmuley and Sohoni tell us to consider partial stability, or (R, P)- stability We must enlarge 
our group action from GL(W), (where W is a m-dimensional subspace of V), to get a larger 
stabilizer of the form: 

/ Mat nxn C * 
^ * 

If we take the natural parabolic, P C GL(W), we have a decomposition, P = KU, where U is 
a unipotent radical and K is the Levi factor. Note that GL(W) is a reductive group. We choose 
a reductive subgroup R C K which has the same rank (in the sense of the rank of the maximal 
torus). Then we choose v € FV such that the following properties hold: 

1) v is instable (The orbit in affine space is closed under the SL(V) action.) 

2) U C G v C P 

Then Mulmuley and Sohoni have a result which is a partial understanding of C[G.v] in terms 
of C[R.v}. Because of the ^-stability, they can gain some understanding from R.v. This stabil- 
ity allows for a significant reduction: Instead of studying the G = GL(n 2 )-module structure of 
C[G.l n ~ m perm m ], we can study the R = SL(m 2 ) -module structure of C[R.perm m ] (See Theorem 
6.2 [44] and the remarks following). 

The addition of the linear factor l n ~ m is also related to the Weyman-Kempf geometric technique 
of calculating syzygies, [58]. To this end one notices that the subvariety 

Y = {/ € S d (W) \f = I s g} 

where I is a linear form, can be treated by means of the geometric technique. One just has to 
consider the incidence variety 

Z = {(f,l)eS d (W)xF(W) | f = l s g}, 

which gives a desingularization of the orbit closure. Here we identify the projective space P(W) 
with the set of non-zero linear forms in W (modulo the usual equivalence relation). There are 
technical difficulties, as Y is usually not projectively normal. 

4. Open Problems 

(1) In the case n = 3, compute the difference between Gl n 2.det n and Mat n 2 xn 2.det n . One 
containment is obvious. In order to understand why the other containment could fail, 
consider the following example: The set Mat2x2-( xd + y d ) does not contain the GL2-orbit, 
GL2-{x d + y d ) for d > 4. This example shows that one must be careful with orbits and 
closures as the behavior is not always as expected. 

In general, compare GL n 2.det n to Mat n 2 xn 2.det n . Note that while Gl n 2 is a group, 
Mat„2 xn 2 is an algebraic monoid. 

(2) (Landsberg) If the two orbits in (1) are not equal, is there a polynomial p(n) such that 
Mat p ( ri )2 xp ( n )2.detp( n ) contains l p( ^~ n Gl n 2 .det n for some linear form 11 

The affirmative answer to this question would imply that VP = VP. 

(3) Let R = Sym(A <g> B <g> C) = A;/1)i , (S\A S^B <g> S v C) Kx ^ , with K x ^ the Kronecker 
Coefficients. Let A(A, B, C) = {(A, fi, v) \ 3n such that S n \A <S> S n)J ,B <8) S nu C C R}. It is 
known that A(A, B, C) is a rational polyhedral cone. Determine the facets of this cone. 

Some cases of this question are already known. For instance in [28], Franz considered 
the case when dim(A) = dim(B) = dim(C), and completely determined the facets in the 
3x3x3 case. His method does not necessarily generalize to give all the facets of the cone. 
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There are also some partial results by Manivel [39], Klyachko [33], and Bernstein-Sjamaar 
[7]- 

(4) The case when dim(A) = 3 and dim(B) = dim(C) = n. Compute the ring of invariants, 

Sym(A®B®C) SL W xSL M 

Note that if dim(A) = 2, then this is the known Sylvester formula, which is related to 
the Kronecker coefficients. This problem would determine K\^n^n for A having < 3 parts, 
(see [40]). 

(5) Consider S\A, where |A| = Yl \ 1S divisible by dim(A), (i.e. A is in the root lattice). Let 
T be the maximal torus of SL(A), i.e. the diagonal matrices with det = 1. Let W denote 
the Weyl group of SL(A). Let (S\A) Q be the zero weight space. Then (S\A) a = (S\A) T 
is a representation of W. For which A does there exist a W invariant? In other words, for 
which A is (S X A)™ + 0? 

Manivel explained that the answer to this question is known to be equivalent to the fact 
that S\A has non zero multiplicity in some S d (S n A). 

Remark 4.1. Questions 3 and 4 relate to the coordinate ring C[GL n 2.det n ], while questions 3 and 
5 relate to the coordinate ring C[G L m 2perm m ]. Question 3 leads to an asymptotic understanding, 
but it is a completely general question. 

Remark 4.2. It is incredible how many important areas come to play. 
Part 3. Matchgates and Holographic Algorithms 

5. REPORT 

In this group, the goal was to understand the work of J.-Y. Cai- P. Lu (see [11] for example) and 
L. Valiant (see [57]) in light of algebraic geometry and representation theory. 

Let F be a field. If T is a planar graph with 7 nodes and F-weighted edges, the FKT (Fisher- 
Kasteleyn-Temperley) algorithm counts the number of perfect matchings for T in polynomial time 
as follows. The algorithm assigns an orientation to the graph (by giving a ±1 to the weight of 
each edge) and associates to T a skew symmetric matrix, k(T). The Pfaffian of k(T) counts the 
number of perfect matchings of T. The Pfaffian and the determinant are known to be computable 
in polynomial time via an application of linear algebra. 

The goal of matchgates and holographic algorithms is to exploit this polynomial time algorithm 
to solve problems which were once unknown to be in P. The methods work when they can change 
a counting problem to one that just requires computing a determinant or Pfaffian. 

An example of the use of matchgates is the following. Consider the problem #7P1-NAE-ICE 
which asks the number (modulo 7) of orientations of a planar graph so that there are no sources 
and no sinks. The first step is to replace all of the vertices and edges of this graph with signatures 
(called generators G and recognizers R) and then try to model these signatures with matchgates. At 
each node, the signature (0, 1, 1, 1, 1, 1, 1, 0) is the truth table values required for the "not all equal" 
condition. At each edge, the signature (0,1,1,0) is truth table values required by the fact that 
each edge carries precisely one orientation. If these signatures can be realized by matchgates, then 
we will be able to compute #7P1-NAE-ICE in polynomial time because this procedure changes the 
question about the old graph to a question about the new graph which can be done in polynomial 
time and gives the desired answer for the old graph. 

For a given problem, what is desired is to know whether there exists a skew symmetric matrix 
associated to a particular signature. This can be determined by examining the so-called matchgate 
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identities (MGI). If the desired signature satisfies MGI, then it is possible to compute that signature 
as the number of perfect matchings of a graph. If the desired signature does not satisfy MGI, all 
is not lost. Special basis changes are allowed as long as they are done in a way that is consistent 
with both the generator and the recognizer - more on this in a moment. 

If / C [7], let Pf 2 i denote the vector (in C 2 ' 7 ') of principal sub-Pfaffians of k(T) produced by 
omitting nodes indexed by J C / and the adjacent edges. G is a standard signature if it occurs as 
a vector of Pfaffians Pf 2 i for some graph T (i.e. the MGI are satisfied). A signature G is realizable 
if there exists a linear map / : C c — > C 2& so that G = f.G. Note that this is just a linear map, and 
not even a group action. The map / is defined by taking a sufficient tensor power of a vector and 
implementing it as a matrix. Similarly, a recognizer signature R is realizable if R = R.f for some 
standard signature R and the same / as before. The goal is to find an / so that both R and G are 
realizable. 

This problem has interesting connections to geometry that the participants hope to be able to 
exploit upon gaining further understanding of these connections. Consider Spec((C 2 *)) k /MGI) = 
$fc. The MGI are well known to geometers, and the varieties that they cut out have been studied 
already. There is a canonical identification of the MGI variety, with the spinor varieties, 
given by the inclusion 

We also have maps / : C c -> C 2 locally defined on the wires of the graph. So, given an explicit 
piece of data, C 2 ®<£? , (assume k < I) and pairing information : [k] <—* [I], and a map / : C c — > C 2 
as above, if c = 2 then the basis collapse theorem of Cai-Lu tells us that we can reduce the k to 2. 
It is unknown what happens if c 7^ 2. 

6. Open Problems 

Exercise: Understand the space of signatures in the product of spinor varieties S^x <S;. What is 
the compatible action? 

(1) Using group actions, reprove the Cai-Lu Basis Collapse Theorem. 

The vector space C 2 <S> C 2 has (at least) two useful descriptions. We can consider the 
tensor product of spin representations, ® A; ~ C 2 * (g) C 2 \ where A p = f\ even C p+1 : this 
reflects the action of the Spin groups on the space. We can also consider the description 
as tensor powers of C 2 , (C 2 )® k <g> (C 2 ) & ~ C 2 " <g> C 2 ': this reflects the action of the 2 x 2 
matrices. 

One may ask: Why not consider a larger space of actions, but the Cai-Lu result reduces 
to studying the action (actually just a mapping) of Mat2(C). (Gurvits comments that the 
spirit of this problem is like quantum computing and unitary operations.) 

The goal of this problem is that by rephrasing the proof in terms of geometry and rep- 
resentation theory, this will make the proof of the basis collapse theorem more transparent 
to geometers, and perhaps it will lead to a greater understanding. It is hoped that this 
greater understanding will lead to a proof of the next problem. 

(2) Study matchgates with questions about colors instead of just yes/no questions. Use the 
results of the previous problem to try to find a natural class of transformation for this 
setup. The idea is that one may be able to use a different group (GL4 for example) in place 
of GL 2 . Does the Cai-Lu Theorem generalize? The expectation is that it will not, and in 
that case, this would be considered good news. This is related to a paper of Valiant which 
considered 3 colors and saw some success. 
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More specifically, if G is a standard generator and R is a standard recognizer, Cai-Lu 
considered the mapping by T, a 2 k x 2 matrix which sends the standard signatures G, R 
to desirable non-standard signatures G, R. What if we consider the case that T is a 2 k x 3 
matrix, or a 2 k x 4, etc, matrix? This would have interesting applications to theoretical 
computer science. 

(3) Consider {£?)® k ®{(C 2 )*)® 1 ~ A^tgiAj. We want to understand the diagonal GL 2 -action on 
points of Sk as an orbit structure. Many cases are already known. For instance, Cai pointed 
out that there is a 2 dimensional space of matrices which leave invariant. What is desired 
is to understand and classify the points which have a 1 dimensional space of matrices which 
move the point off of the spinor variety. The reason is that if a point has this feature, there 
will then be 1 dimension's worth of freedom for the each factor in Seg{Sk x Si). This in 
turn will tell which generators have a 1-parameter family of compatible recognizers. 

Part 4. Entanglement and Quantum Information Theory 

7. Report 

The contributions of G. Gour and D. Gross mention the long-standing additivity conjecture 
of quantum information theory 0. D. Gross posed a question which is closely related to this 
open problem - but has the advantage of being independent of the measure used to quantify 
entanglement. In particular, there is no need to talk about non-algebraic quantities such as the 
"entropy of entanglement" . 

The basic setting of the additivity conjecture is a linear space with a two-fold (horizontal and 
vertical) tensor product structure 

(A! (g>B!) 
® (A 2 ® B 2 ) . 

Physically, the spaces A = A\ ® A 2 and B = B\ ® B 2 are associated with distinct observers, in 
a way that need not concern us. Let S\ be any subspace of A\ (g> B\ and likewise 5*2 C A 2 <S> B 2 . 
Set S = S\ <® S 2 . 

We are interested in the set of vectors in S which display the "least amount of entanglement" 
with respect to the partition A vs. B. The physical quantity "entropy of entanglement" measures 
the deviation of a vector s £ S from being decomposable - of the form sa ® sb- (The precise 
definition will not be relevant for the question asked below). The additivity conjecture amounts 
to claiming that the set of least entangled vectors always contains one which is decomposable with 
respect to the partition 1 versus 2@ In other words, it is conjectured that if one aims to identify the 
global minimum of the entropy of entanglement in S, one may restrict attention to vectors of the 
form si ® s 2 , Si & Si. It is likely that many researchers in the field would agree to slightly stronger 
and less precise statement that "generically" the global minimum is attained only on decomposable 
vectors. 

D. Gross also asks the following consistency question. The optimization problem we set out to 
solve (find the least entangled vector in S) does not refer to the vertical tensor product structure. 
However, the conjectured statement does. So is the vertical decomposition intrinsic to the problem? 

^Soon after the workshop, the additivity conjecture for minimum output entropy was proven false by M.B. Hastings, 
[29]. 

The conjecture pertains only to the case where the entropy of entanglement is used to measure the degree of 
non-decomposibility of a vector. Analogous statements for several other measures - most notably the tensor rank - 
have recently been proven to be wrong. 
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Given just the data A, B and S C A®B, when can one introduce such an additional tensor product 
structure and when is it unique? 

If a negative answer could be given - i.e. if the tensor product structure is essentially uniquely 
determined by S - we would have learned an interesting fact about tensor product spaces. A 
positive answer would be an extremely interesting starting point for the search of a counterexample 
to the additivity conjecture. 

Deciding additivity is considered the most high-profile open problem of quantum information. 
Any partial result would be met with great interest. There is no compelling evidence either in favor 
or against the conjecture. 

The group working on entanglement and quantum information theory reports that while they 
did many calculations, they did not solve any of the open problems. 

In addition, they computed for s £ S181S2, the quantity mm{H(s)}, where H can be the entropy 
or the non-locality for example. They claim that it is also possible to compute 



and in so doing, they find s° p <8> s^f . It is unknown whether this is a local optimum. The following 
statement has been solved: If s° pt <g> is a local optimum, s opt ® s^ 1 is a critical point, and this 
implies that locally to first order that additivity works. Unfortunately, this tells nothing about the 
global problem. 



(1) To which degree can one violate the Rank Multiplicitivity (simpler version)? 
Note: This problem is perhaps not so interesting as it was aimed at the additivity conjecture for 
minimum output entropy which was recently proved to be false. Nonetheless, we report on the 
activity from the workshop surrounding this problem. 

Consider subspaces S\ C A\ ® A2, S2 C B\ ® B2, and let r(Si) = min ^ sg 5 i {ranA;(s)}. Let 
S = S\ <g> S2. We may consider S C {A\ <g) A2) <8> {B\ <g> B2), or under the natural braiding 
isomorphism, we can consider S C (A\ <8> B{) <g> {A2 <8> B2). It is known that there exist examples 
where r(S) 7^ r(5i)r(S , 2). Find examples of large violations. (Make it precise what it means to be 
"large".) It should be noted that this is a toy example. There is not an automatic correspondence 
to the big problem for entropy or entanglement. 

Following are two examples which were worked out during the workshop. Let X C Mat nxn (R) 
be a linear subspace of real valued matrices, and let Y = X®X. Define MinRank(X) as a minimal 
rank of nonzero matrices in X. Clearly (MinRank(X)) 2 > MinRank(Y). 

Proposition 8.1 (L. Gurvits). The MinRank - multiplicativity conjecture, (MinRank(X)) 2 = 
MinRankiY) , is false overR. Moreover, there can be an arbitrarily large decrement (MinRank(X)) 2 — 
MinRank(Y)) overR. 



X = span{M,l2}. Clearly, the minimum rank of any nonzero element in X is 2. The spectra 
a{M®M) = (1,1,-1,-1) and M®M is symmetric. Therefore, we find that rank{M M + I A ) = 
rank(M (g> M — I 4) = 2 < 4. So MinRank multiplicativity is false. 

For the arbitrary large decrement, modify the example as follows. Construct the linear space 
S2n = X ® I n = span{M <8> I n , I2 <2> In = l2n}- Then the minimal rank of nonzero matrices in 
S2n is 2n, but rank((M <g) I n ) <g) (M <8> I n ) — I 4^) = In 2 and this provides an arbitrarily large 
decrement. □ 



S16S1 s 2 £S 2 



min H(s±) + min H(s2), 



8. Open Problems 
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Remark 8.2. S. Friedland observed that the same example shows that the entropy additivity con- 
jecture is false over reals since log(2n) + log(2n) > log(2n 2 ). 

(2) Let dim(Si) = fcj. Assume each Si is generic in the sense that it is in general position in 
Ai (g) Bi. Is it the case that you can have r(S) = r(Si)r(S2)7 This is not just a dimension 
counting argument because S± £3 S2 may not be generic. In fact, it was commented that 
S\ <8) S2 is never generic. (See Gurvits's example for the real case above). 

(3) Ambiguity of Tensor Decomposition. This area hasn't been explored much. 

Are there interesting situations which allow us to find a tensor product structure A = 
A[®A' 2 , B = B[®B' 2 , (A' i7 t A u B\ / Bi) and subspaces S' t C A^®B' t such that S = S[<g>S 2 ? 
To qualify as "interesting", S should contain no unentangled vector (of the form sa <8> s_b). 
Also, it would be desirable if the two decompositions A = A\ ® A2 and A = A'-y ® A 2 shared 
no decomposable vector. (There should be some more precise conditions given for a solution 
to be interesting). 

The uniqueness of this decomposition is interesting. If you prove non- uniqueness, then 
you become a celebrity as it is related to an important problem in entanglement and quan- 
tum information theory. 
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